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ABSTRACT 

We present evidence from a sample of 544 galaxies from the DEEP2 Survey for evolution of the inter- 
nal kinematics of blue galaxies over 0.2 < z < 1.2. DEEP2 provides a large sample of high resolution 
galaxy spectra and dual-band Hubble imaging from which we measure emission-line kinematics and 
galaxy inclinations, respectively. Our large sample allows us to overcome scatter intrinsic to galaxy 
properties, in order to examine trends. At a fixed stellar mass, galaxies systematically decrease in 
disturbed motions and increase in rotation velocity and potential well depth with time. The most 
massive galaxies are the most well-ordered at all times, with higher rotation velocities and less dis- 
turbed motions compared to less massive galaxies. We quantify disturbed motions with an integrated 
gas velocity dispersion (cr 9 ), which is unlike the typical pressure-supported velocity dispersion mea- 
sured for early type galaxies and galaxy bulges. Due to finite slit width and seeing, a g integrates over 
unresolved velocity gradients which can correspond to non-ordered gas kinematics such as small-scale 
velocity gradients, gas motions due to star-formation, or super-imposed clumps along the line-of-sight. 

We compile surveys of galaxy kinematics over 1.2 < z < 3.8 and do not find any trends with redshift, 
likely because these studies are biased toward the most highly star-forming systems. In summary, 
over the last ~ 8 billion years since z = 1.2, blue galaxies evolve from disturbed to ordered systems as 
they settle to become the rotation-dominated disk galaxies observed in the Universe today, with the 
most massive galaxies always being the most evolved at any time. 

Subject headings: galaxies - formation, galaxies - evolution, galaxies — kinematics and dynamics, 
galaxies - fundamental properties 


1. INTRODUCTION 

In the standard picture of disk galaxy formation, 
baryons enter into dark matter halos and dissipate to 
form disks at their centers (Rees & Ostriker 1977; White 
& Rees 1978; Fall & Efstathiou 1980; Blumenthal et al. 
1984). The accretion of baryons onto galaxies is expected 
to be an on-going process, and it is not yet known how 
galaxy disks respond to it, or how they evolve as mass ac- 
cretion rates change. Only recently have hydrodynamic 
simulations begun to address this (Kimm et al. 2011; 
Danovich et al. 2011). The cleanest observational probe 
of disk galaxy assembly is internal galaxy kinematics. It 
tells us directly about the dynamical state of galaxies, re- 
veals the potential well depths of individual galaxy-dark 
matter halo systems, and can be measured for a large 
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sample of galaxies over a significant range in redshift. 

We have acquired a wealth of information about blue 
galaxies over the last ~ 8 billion years since z — 1.2 which 
can be roughly divided into two categories: those which 
imply a significant amount of evolution and those which 
do not. We posit that the best way to resolve these seem- 
ingly contradictory results is with a study of the internal 
kinematics of galaxies. Of those studies which do not 
imply much evolution are galaxy luminosity functions. 
They find that the number density of blue galaxies is 
unchanging, and that blue galaxies fade by only ~ 1 B- 
band magnitude since z ~ 1 (e.g., Willmer et al. 2006; 
Faber et al. 2007; Bell et al. 2005). Similarly, stellar 
mass functions of galaxies do not evolve to within uncer- 
tainties on the stellar mass measurements (~ 0.2 dex), 
although there are some hints for a slight decrease in 
number density with time (Drory et al. 2009). In addi- 
tion, galaxy sizes are only marginally smaller at z = 1 
(by a factor of 1.4; Dutton et al. 2011), and the major 
merger rate of blue galaxies is only ~ 10 - 20% between 
z — 1 and now (Lotz et al. 2011). However, despite all 
this evidence for a slowly evolving population, there are 
some strong indications to the contrary. Namely, blue 
galaxies at z ~ 1 have significantly higher star-formation 
rates (by a factor of ~ 3; e.g., Noeske et al. 2007), higher 
molecular gas fractions for at least the brighter galax- 
ies at z = 1.1 — 1.2 by a factor of ~ 2 - 4** (Daddi et 
al. 2010; Tacconi et al. 2010), and more disturbed mor- 
phologies (e.g., Abraham et al. 1996; Abraham & van den 
Bergh 2001). In summary, galaxies do not evolve much 
in luminosity, stellar mass, or size, but they do evolve 
strongly in star-formation rate, molecular gas mass, and 
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morphology. It is unclear what mechanisms are behind 
these transformations. In this paper, we attempt to un- 
derstand this discrepancy, and the mechanisms behind 
it, by studying the internal kinematics of blue galaxies 
over 0.2 < z < 1.2. 

Due to the large intrinsic scatter in the kinematic prop- 
erties of blue galaxes (Weiner et al. 2006a, b; Kassin et 
al. 2007), a sizeable sample (> 100) which is represen- 
tative in terms of galaxy properties is required to study 
internal galaxy kinematics over a significant look-back 
time. Such a sample is provided by the Deep Extragalac- 
tic Evolutionary Probe 2 Survey (DEEP2; Newman et al. 
2012). DEEP2 is one of two large galaxy surveys which 
cover a significant range in redshift and have high enough 
spectral resolution to measure semi-resolved galaxy kine- 
matics, the other being the much smaller Team Keck 
Redshift Survey (TKRS; Wirth et al. 2004). These sur- 
veys were created with galaxy kinematics in mind. In 
addition, a portion of DEEP2 galaxies also have Hubble 
Space Telescope/ ACS imaging in two bandpasses (V and 
7; F606W and F814W, respectively). Hubble images are 
required for measurements of galaxy axial ratios, which 
are needed to correct rotation velocities for inclinations 
of galaxies to the line-of-sight. The two wavebands al- 
low for axial ratios to be measured at similar rest-frame 
wavelengths for galaxies over 0.2 < z < 1.2. 

The DEEP2 and TKRS Surveys have revealed that the 
kinematics of blue galaxies are not fully described by ro- 
tation velocity (Vi-oi) alone (Weiner et al. 2006a, b; Kassin 
et al. 2007). About ono-third to one-half of emission-line 
galaxies over 0.2 < z < 1.2 have a significant or domi- 
nant component of disordered motions, as measured via 
an integrated gas velocity dispersion ( a g ). The impor- 
tance of cr s in galaxy kinematics was also found by sur- 
veys using Integral Field Spectroscopy (IFS) of represen- 
tative galaxies at z ~ 0.6 (IMAGES; Yang et al. 2008), 
and strongly star-forming galaxies at z ~ 1.3 (Wright et 
al. 2007, 2009; Epinat et al. 2009; Lemoine-Busserolle & 
Lamareille 2010; Vergani et al. 2012), z ~ 2 (Forster- 
Schreiber et al. 2009; Law et al. 2009), and z ~ 3 
(Lemoine-Busserolle et al. 2010; Gnerucci 2011). 

Galaxies with large values of <r g scatter to low V TOt 
in the Tully-Fisher relation (V rot versus luminosity or 
stellar mass) (e.g., Kassin et al. 2007). However, when 
Og is incorporated along with V rot into a kinematic in- 
dicator which is created to trace the pote ntial wells of 
galaxy-dark halo systems, namely Sb .5 = ^/0.5V/ ot + 
a more complete measure of the internal motions of 
galaxies results. When Vrot is replaced by So . 5 in the 
Tully-Fisher relation, there is significantly reduced scat- 
ter (Weiner et al. 2006a; Kassin et al. 2007; Catinella 
et al. 2011). This tightening of Tully-Fisher once S0.5 is 
adopted was also found for IFS observations (Cresci et al. 
2009; Puech et al. 2010; Lemoine-Busserolle & Lamareille 
2010; Lemoine-Busserolle et al. 2010; Vergani et al. 2012), 
and numerical simulations of interacting galaxies (Cov- 
ington et al. 2010). It is clear from these, studies that a g 
captures a necessary component of kinematic support, 
although it is unknown what phenomena cause it. Possi- 
bilities include small-Bcale velocity gradients, gas motions 
due to star-formation, and super-imposed clumps along 
the line-of-sight. 

Kassin et al. (2007) showed that galaxies evolve in V ro t 
and (J g since z = 1.2, and quantified the evolution in 
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Fig. 1. — Rest-frame U — B color versus stellar mass diagrams for 
galaxies with solid spectroscopic redshifts in field 1 of the DEEP2 
Survey (grey points) and for the K07 sub-sample (black points) are 
shown in bins of redshift. The K07 sample traces the full stellar 
mass range of the DEEP2 Survey and a significant range in color 
in the “blue cloud” at all redshifts. Dotted lines demarcate the 
mass-complete sample used later in this paper. 


terms of the Tully-Fisher relation. In this paper, we 
study the evolution of V ro t, o g and S 0.5 individually. 
We ask the following questions of our dataset: Do disks 
form/settle from an initially disordered state? If so, what 
are the kinematic signatures of disk settling/formation? 
Is this process dependent on galaxy properties such as 
stellar mass and size? What is the fraction of blue galax- 
ies which are settled today and at z = 1? A key element 
of this study is our broad sample selection, which is de- 
scribed in §2. In §3 we discuss the stellar mass and size 
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Fig. 2.— Relationships between kinematic quantities (<r 9 , V r „t, Sc. 5) and stellar mass are shown in bins of redshift for the K07 sample. 
While large scatter is found for the a g and V ro t relations, the So . 5 relation is significantly tighter and does not evolve to within uncertainties 
since 2 =' 1.2. The V rot and So.s relations are the same as in K07; we Bhow the a g relation here for completeness. Hashed regions denote 
areas where the survey is not sensitive, and dotted lines demarcate the mass-complete sample used later in this paper. If a galaxy is deemed 
too morphologically disturbed to accurately determine an inclination (9% of the sample), Wot is not inclination corrected, and the galaxy 
is shown as an upper limit symbol in the V rat and S 0.6 plots and an open circle in the a g plot. 


measurements adopted, and in §4 we describe how Vrot 
and a g are measured from emission lines in galaxy spec- 
tra. We show how V rot ,a g , and So . 5 evolve with redshift 
in §5, and discuss the influence of galaxy stellar mass 
and size on this evolution in §6. In §7 we quantitatively 
define a settled disk galaxy and study their fraction with 
redshift. We compare our results to measurements of 
galaxy kinematics in the literature over 0 < z < 3.8 in 
§8. Conclusions are given in §9. A ACDM cosmology is 
adopted throughout (ft = 0.7, fi m = 0.3, 0^ = 0.7), and 
all logarithms are base 10. 

2. SAMPLE SELECTION 

One of the main differences between this paper and 
most other studies of galaxy kinematics lies in how our 
sample is selected. Typically, constraints on galaxy mor- 
phology are applied to remove disturbed and major- 
merger galaxies. We instead include all galaxies with 
emission lines which are bright enough for us to mea- 
sure kinematics from, regardless of morphology. Since 
disturbed morphology and disordered motions (as mea- 
sured by o g ) are correlated, removing such galaxies in- 
troduces a strong bias toward only the most rotation- 
dominated systems (Kassin et al. 2007). Little or no evo- 
lution in internal galaxy kinematics would be found for 


a sample selected in this manner. Such a sample would 
not be representative of galaxies at higher redshifts since 
they tend to have disturbed morphologies (e.g., Abra- 
ham et al. 1996; Abraham & van den Bergh 2001). Even 
when a strict morphological selection is loosened some- 
what to include galaxies which are moderately disturbed 
but disk-like, still little scatter and no evolution in Tully- 
Fisher are found (Miller et al. 2011, 2012). Furthermore, 
many studies of galaxy kinematics are limited to galaxies 
with extremely bright emission lines. The sensitivity of 
DEEP2 allows us to probe a large range in galaxy stellar 
mass and color (Figure 1). 

In addition to a broad morphological selection crite- 
rion and a large range in stellar mass, in order to study 
the evolution of disk galaxy kinematics over a significant 
look-back time, high resolution spectra of a sizable sam- 
ple of galaxies which span a wide range in redshift are re- 
quired. A sizable sample is needed because emission-line 
galaxies have significant intrinsic scatter in their kine- 
matic properties (Weiner et al. 2006a, b; Kassin et al. 
2007). Tb this end, we adopt measurements of inter- 
nal galaxy kinematics for 544 galaxies from Kassin et al. 
(2007) (hereafter K07). This sample derives from a very 
large collection of 9715 galaxies with solid spectroscopic 
redshifts from Keck/DEIMOS (Faber et al. 2003) in field 
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Fig. 3. — Left: The ratio of ordered to disordered motions in a galaxy, V Tot /<x s , is shown as a function of stellar mass and in bins of 
recshift for the K07 sample. The galaxies define an upper envelope which does not evolve with redshift. However, galaxies evolve within 
the envelope: the lower-right corner evacuates with time (from high to iow redshift), with higher mass galaxies evolving the fastest. Hashed 
regions denote areas where the survey is not sensitive, and dotted lines demarcate the mass-complete sample used later in thiB paper. 
Galaxy symbols are as in Figure 1. Right: This evolution is more clearly shown as histograms of Vrat/og for galaxies divided into stellar 
mass bins. 


1 of the DEEP2 Survey (May 2007 catalog). Galaxies We now detail how galaxies were selected from field 1 
in field 1 of DEEP2 are observed to a limiting magni- of the DEEP2 Survey to be in the K07 sample. First, 

tude of Rab <24.1 and axe not subject to a color se- although the DEEP2 reaches a redshift of z = 1.5, K07 

lection (Willmer et al. 2006). In Figure 1, we show rest- chose to be conservative in their selection and limit their 

frame color-magnitude diagrams for galaxies in field 1 of sample to z = 1.2, leaving 8787 galaxies. This helps to 

DEEP2 and the K07 sample. The K07 sample traces the avoid the high-redshift end of the survey where the se- 

entire stellar mass range of the survey and a significant lection is in the rest-ultraviolet and dependent on color 

range in color in the “blue cloud” at all redshifts. (Willmer et al. 2006). It also avoids galaxy axes ratios 
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measured from Hubble imaging which has shifted into 
the ultraviolet (axes ratios are used to correct rotation 
velocities for inclinations of galaxies to the line-of-sight). 
Next, galaxies were chosen to have Hubble Space Tele- 
scope/ ACS imaging at V and J, which further reduces 
the sample to 3523 galaxies. This is because there are 
portions of field 1 of the DEEP2 Survey which are not 
imaged by Hubble. Beside removing galaxies which did 
not have Hubble imaging, this requirement did not affect 
our sample selection any further because the Hubble im- 
ages are deeper than the spectral survey. For our most 
important selection, galaxies with bright enough emis- 
sion lines to measure kinematics (> 10“ 17 erg s' 1 cm 2 ) 
were chosen, leaving 1692 galaxies. 

Galaxies were further selected to have inclinations 
measured from Hubble images between 30° < i < 70°, 
and spectrographic slits with position angles aligned to 
within 40° of their major axes. These cuts leave a sample 
of 755 galaxies. Inclinations are measured from Hubble 
images using the GIM2D software (Simard et al. 2002), 
as described in Dutton et al. (2011) and Cheung et al. 
in preparation. The inclination cut avoids nearly face- 
on galaxies for which inclination measurements are very 
uncertain, and the effects of dust in the determination 
of stellar mass for highly inclined systems. Inclinations 
are measured from the V-band Hubble image for galax- 
ies at z < 0.6, and the /-band Hubble image for galaxies 
at z > 0.6. However, if the morphology of a galaxy 
was deemed by eye to be disturbed enough such that its 
inclination could not be reliably determined, it was in- 
cluded in the sample but not corrected for inclination. 
There are only 24 of these galaxies and they are flagged 
in the following analysis. Removing them from the sam- 
ple does not significantly affect the results of this paper. 
For the slit position angles, due to the large slit width 
(1") compared to the apparent sizes of galaxies in our 
sample (~ 3"), seeing will allow us to measure accurate 
kinematics when slits have position angles which are off- 
set by up to 40° from galaxy major axes (see Figure 13 of 
Weiner et al. 2006a). Finally, if an emission line was af- 
fected by a sky line or instrumental artifact, we removed 
the corresponding galaxy from the sample, leaving a to- 
tal of 544 galaxies in the K07 sample. No constraints on 
morphology were ever applied. 

To study the evolution of galaxy kinematics for a sam- 
ple which is mass-complete at all redshifts, in §5.1 we 
will limit the K07 sample to galaxies with stellar masses 
(A/.) over 9.8 < logM.(MQ) < 10.7, as demarcated in 
Figures 1- 3. Since the M, function of blue galaxies does 
not evolve significantly since z ~ 1 (Borch et al. 2006; 
Bundy et al. 2006; Pozzetti et al. 2010), we do not change 
this M, range with redshift. This cut removes 274 galax- 
ies from the K07 sample, leaving 270. Throughout the 
paper we will refer to this as the “mass-complete sam- 
ple” 

3. STELLAR MASS AND SIZE MEASUREMENTS 

In this paper, we study the relationship between galaxy 
kinematics and M.. We adopt M, measurements from 
Bundy et al. (2006) who performed fits of galaxy stel- 
lar population models to observed B, R, I, and iZ„-band 
integrated galaxy fluxes. The models constituted a grid 
of synthetic spectral energy distributions from Bruzual 
& Chariot (2003) which spanned a range of exponential 
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Fig. 4. — For the mass-complete sample, the distribution of 
galaxy stellar masses among the redshift bins in Figure 3 is shown. 
It demonstrates a bias in the lower redshift bins towards lower 
stellar masses, and a bias in the higher redshift bins towards higher 
stellar masses. These biases work in the opposite direction to our 
results, making our results lower limits to the intrinsic evolution. 


star-formation histories, ages (restricted to be less than 
the age of the Universe at the redshifts of the galaxies), 
metallicities, and dust contents. A Chabrier (2003) ini- 
tial mass function was adopted. For a given galaxy, at 
each point on the grid of models, the following was cal- 
culated: the minimum y 2 , probability that each model 
accurately describes the galaxy, and the /fa-band stellar 
mass-to- light ratio. The corresponding M, at each point 
was then determined by scaling the stellar mass-to- light 
ratio to the AVband luminosily. The probabilities were 
then summed (marginalised) across the grid and binned 
by model M *, which yielded a M, probability distribu- 
tion for the galaxy. The median of the distribution was 
adopted as the best estimate. The M„’s measured in 
this way are robust to degeneracies in the models, such 
as those between age and metallicity. The uncertainties 
in these M m measurements are ~ 0.2 dex. 

In addition to the effects of M, on galaxy kinemat- 
ics, we also look to galaxy sizes. Two size measurements 
are studied: emission line extents and continuum sizes. 
Emission line extents are measured from the ground- 
based DEEP2 spectra from the same lines used to mea- 
sure kinematics. However, unlike the kinematics mea- 
surements, seeing is not corrected for. We measure the 
spatial extent of the line emission as follows. For de- 
tails on the method, see Weiner et al. (2006a). First, 
we measure the spectral continuum by collapsing the 
two-dimensional spectra in wavelength over a 15 A range 
around the line and two 100A spatial ranges on either 
side of the line. We then fit Gaussian profiles along the 
slit using a non-linear least-squares routine, and tabulate 
the resulting full-width-half-maxima (FWHM) , which we 
adopt as our measurement of emission line extent. Our 
other size measurements, galaxy continuum sizes, are 
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measured from Hubble images, and are therefore unaf- 
fected by seeing, lb take into account band-shifting with 
redshift, sizes are measured from the V and /-band im- 
ages for galaxies at redshifts less than and greater than 
1.6, respectively. The images are of one-orbit depth and 
reach 28.75 and 28.10 AB magnitudes (5<r) for V and /, 
respectively. Measurements of galaxy half-light radii are 
made with the GIM2D software (Simard et al. 2002), as 
described in Dutton et al. (2011) and Cheung et al in 
preparation. 

4. KINEMATIC MEASUREMENTS 

In this section, we describe the measurements of in- 
ternal galaxy kinematics from emission lines which we 
adopt from K07. Details on how these measurements 
were made are in Weiner et al. (2006a) for a similar sur- 
vey with the same telescope/instrument, but at lower 
spectral resolution ( R ~ 2100 versus R ~ 5000). The 
emission line used to measure kinematics is the one with 
the highest signal-to-noise in the galaxy spectrum. For 
the vast majority of galaxies, the Ha A6563, [O II] 
AA3726.0, 3728.8, and [O III] A5007 lines are used. Ro- 
tation velocities on the flat part of the rotation curves 
(V ro t xsm(i), where i is the inclination of a galaxy to the 
line-of-sight) and integrated gas velocity dispersions (o g ) 
were measured simultaneously from emission lines in the 
spectra, and the effects of seeing (typically ~ 0.7") were 
taken into account. Our high spectral resolution allows 
for V rot and o g to be measured down to ~ 5 and ~ 10 
km s -1 , respectively, if present. An uncertainty of 10 km 
s -1 is adopted for both a g and the inclination-corrected 
V rot to account for random errors and the dependence of 
the model parameters on the assumed seeing and scale 
radius of the rotation curve. Although many of the rota- 
tion curves do not “turn over” due to seeing, our model 
is still able to fit for V rot . It has been demonstrated us- 
ing data with the same telescope/instrument which reach 
the turn-over radius for 90% of galaxies, that spectra of 
our depth do not show a bias in V rat (Figure 6 of Miller 
et al. 2011). However, benefits of the deeper data are 
smaller errors in V ro t x sin(i) (and therefore less scatter) 
and the ability to probe galaxies with fainter emission 
lines (Miller et al. 2011). Except for the 24 galaxies for 
which reliable inclinations could not be determined (§2), 
values of V^t x sin(i) are corrected for galaxy inclina- 
tions measured from the Hubble images. The results of 
this paper do not change significantly if these 24 galaxies 
are removed from the sample. 

4.1. What o g and V ra t Measure 

The integrated gas velocity dispersion we measure (o g ) 
is unlike the typical pressure-supported velocity dispersion 
measured for early type galaxies or galaxy bulges. This is 
the case for two reasons. First, o g is measured from emis- 
sion lines which trace the gas in galaxies, as opposed to 
absorption lines which trace stars. Gas, unlike the colli- 
sionless stars in an early type galaxy, can dissipate energy 
and therefore cannot remain in a high dispersion equilib- 
rium state with crossing orbits. Secondly, because the 
spectral slits used in the DEEP2 Survey axe wide com- 
pared to the apparent sizes of the galaxies observed (1" 
slit vs. ~ 3" galaxies), o g effectively integrates velocity 
gradients on scales at and below the seeing limit (Weiner 
et al. 2006a; Covington et al. 2010; Epinat et al. 2010). 


For values of o g less than ~ 35 km s -1 , which is the upper 
limit for well-ordered disk galaxies in the local Universe 
(Figure 15 of Epinat et al. 2010), o g measures the relar 
tive motions of individual star-forming regions in galaxy 
disks. However, for values of o g greater than ~ 35 km 
s -1 ,- o g integrates over velocity gradients which can corre- 
spond to non-ordered gas kinematics such as small-scale 
velocity gradients, gas motions due to star-formation, or 
superimposed clumps along the line of sight. In this pa- 
per, we collectively refer to these velocity gradients as 
“disturbed motions. ” 

To test for interference from V ra t in measurements of 
o g , in §5 we compare o g measured for galaxies in our 
mass-complete sample with galaxies from DEEP2 which 
are face-on in Hubble images (i.e., i < 30°) and in the 
same M« range. We find no significant difference in the 
trend of o g with redshift between these two samples, 
making it unlikely that the V ro t of a galaxy interferes 
with the measurement of o g . 

To further investigate what o g and V rot measure, in 
Covington et al. (2010) we performed mock observations 
of a suite of gaiaxy merger simulations. These simula- 
tions provided model galaxies in various stages of disor- 
der, ranging the gamut from well-ordered disks to merger 
remnants. The mock observations were performed just 
as the actual observations. Galaxies were redshifted to 
z ~ 0.3 and z ~ 1.0, and seeing, slit width, and detector 
pixel size were taken into account. Our algorithm for 
fitting kinematics was found by Covington et al. (2010) 
to be successful at reproducing o g , but was found to un- 
derestimate values of V Tat by up to 30%, independent of 
a merger stage or redshift. However, as we note in §2, 
Miller et al. (2011) compared measurements of V TOt for 
data similar to ours with deeper data which reach the 
turn-over radius, and do not find a systematic offset in 
V ro t between the two measurements. Therefore, if we 
are underestimating V ro t, we are probably not doing so 
significantly. 

5. EVOLUTION OF GALAXY KINEMATICS 

In this section, we examine how the kinematics of blue 
galaxies evolve with redshift. First we examine the full 
K07 sample, and then we study a mass-complete sub- 
sample. 

5.1. The Full K07 Sample 

In Figure 2, we show o g , V r ^, and So.s versus M. in 
bins of redshift for the K07 sample. The V rat versus 
M, relation (i.e., the stellar mass Tully-Fisher Relation) 
shows large scatter to low V TO t (K07). The o g versus M» 
relation also shows significant scatter in o g . However, 
when V^ and o g are combined into S0.5, the resulting re- 
lation with M« tightens significantly and does not evolve 
with redshift to within uncertainties (K07). The quan- 
tity S0.5 is a tracer of potential well depth (Weiner et 
al. 2006a; Kassin et al. 2007). Figure 2 also shows that 
at lower redshifts, higher mass galaxies have on average 
larger values of V ro t and smaller values of o g than at 
higher redshift. It is also clear from this figure that at 
lower redshift, lower mass galaxies have low values of V ro t 
and high values of o g ; they are too faint to be observed 
at higher redshifts in our survey. 

These findings are illustrated more clearly in Figure 3 
(left column) where the quantity V ro t/x g . the ratio of 
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Fig. 6. — The evolution of a 3 , V ro t, and Sa.s = ^/O.SV^ + with redshift for the mass-complete sample (9.8 < log M.(Mq) < 10.7) is 
shown. Galaxies increase in tr g , and decrease in V ro t arid So . 5 with increasing redshift. Galaxies are plotted as black symbols, as described 
in Figure 1. Binned medians are shown as red triangles with error bars representing the error on the medians determined by bootstrap 
re-sampling the data. Linear fits to the binned medians are shown as red solid lines. 



Fig. 6. — The evolution of Oj/So.s, V'rot/So. 5 , and V ra t/oji with redshift for the mass-complete sample is shown. The first two of these 
ratios remove the dependence on potential well depth (i.e., So s) to the a g and Vrot plots in Figure 5. Symbols, error bars, and fits are the 
seme as in Figure 5. 


ordered to disordered motions in a galaxy, is plotted ver- 
sus M m in bins of redshift. 12 The galaxies define an up- 
per envelope in V ra t/o g versus M.-space which does not 
evolve since z — 1.2. However, galaxies evolve within 
the envelope: with increasing time they move away from 
the bottom of the plot toward the upper envelope, and 
this happens first for massive galaxies. This evolution 
is better quantified by the histograms in Figure 3 (right 
column) which demonstrate that galaxies move towards 
higher V rot l<J g decreasing redshift (i.e., with time), and 
the highest mass galaxies have on average the highest 
values of V ra t/(T a at all times. 

5.2. A Mass-Complete Galaxy Sample 

lb investigate this trend oiV Tot ja g with redshift for all 
but the very lowest and very highest mass galaxies in our 
sample, we create a mass-complete sample by selecting 
galaxies with stellar masses over 9.8 < log M,(M®) < 
10,7 (vertical dotted lines in Figures 1-3). For this sub- 

15 In the evaluation of V Ta t/a 3 , a minimum value of <r 9 of 15 km 
s" 1 is adopted to avoid dividing by small numbers with relatively 
large uncertainties. 


section and for §6, we will focus on this mass-complete 
sample of 270 galaxies, unless noted. In Figure 4, for 
the mass-complete sample, we show the distribution of 
galaxy M, among the redshift bins in Figure 3. It il- 
lustrates that the lower redshift bins are biased towards 
lower mass galaxies, and the higher redshift bins towards 
higher mass galaxies. These biases work in the opposite 
direction to our results. Therefore, the evolution we find 
is interpreted as a lower limit to the intrinsic evolution. 

Figure o shows the evolution of cr 9 , V rot and So.s with 
redshift for the mass-complete sample. All three quan- 
tities have large intrinsic scatter. We will demonstrate 
that a portion of the scatter for V Ta t and o g is due to M., 
or size, in §6. We fit linear relations to each of the quan- 
tities versus redshift by performing linear least-squares 
fits to the binned medians in Figure 5. The medians are 
zero-pointed near the middle of the sample to avoid co- 
variances. Of the three quantities, o g shows the strongest 
evolution (5.0cr). It increases linearly with redshift to 
z = 1.2 (i.e., decreases with time) as 

log <r 9 — 0.42 = (0.25±0.05)(z— 0.82)+(1.16±0.01). (1) 
The median increases from 27 ± 1 km s -1 at z = 0.2 to 









8 


Kassin et al. 



Fig. 7. — The evolution of a 9 with redshift for galaxies oriented 
face-on in Hubble images (» < 30°) is compared with that for 
galaxies in the mass- complete sample from Figure 5, and is not 
found to be significantly different. The solid line is a fit to the 
face-on galaxies, and the dashed line is the fit to the mass-complet e 
sample. Both samples span the same M, range, and both fits are 
given in the text. Symbols and error bars are as in Figure 5. 

47 ± 1 km s -1 at z — 1.2. The quantity which shows 
the second-most strongest evolution is V rat (4.2«x). It 
decreases with redshift as 

log Kot -2.00 = (-0.21 ±0.05)(z — 0.82) — (0.04 ±0.02). 

, (2) 

The median decreases from 123 ± 1 km s 1 at z = 0.2 to 
76 ±1 km s _1 at z = 1.2. Finally, the quantity S0.5 shows 
the weakest evolution of the three (3.6<r). It decreases 
with redshift as 

log S 0 5 - 1-92 = (-0.11 ±0.04)(z - 0.82) + (0.00 ± 0.01). 

1 (3) 

The median decreases from 97 ± 1 km s -1 at z = 0.2 to 
75 ± 1 km s -1 at z = 1.2. 

Some of the evolution found for these three quantities, 
may he due to a changing distribution of galaxy masses 
in the sample with time. This can be caused by galax- 
ies growing in mass and/or high mass galaxies quenching 
their star-formation and leaving the sample. To address 
this for <j g and V rot , we remove the dependence on po- 
tential well depth by examining ratios with So. 5, namely 
Og /So.s and Vrat/So.s- These ratios are shown as a func- 
tion of redshift in Figure 6. They express the fraction 
of kinematic support which comes from ordered or disor- 
dered motions. The ratio ff 9 /So.s increases with redshift 
(5 .0a) as 

log ag/So.o +0.35 = (0.45±0.09)(z-0.82) + (0.01±0.03), 

( 4 ) 

and the ratio V rot /So.s decreases with redshift (3.0<r) as 
log V rot /S 0 5 -0.10 = (—0.06 ±0.02)(z—0.82)+(0..00 ±0.01) 

( 5 ) 

As before, these are linear least-squares fits to the median 
points in Figure 6. These ratios demonstrate that the 


increasing role of <r g and the declining role of V rat with 
increasing redshift are not due to a changing distribution 
of galaxy masses in the sample with time. 

In addition to ratios with S0.5, we look to V TO t/o g . the 
ratio of ordered to disordered motions, with redshift 12 . 
This is not an ideal quantity since neither V ra t nor er g is 
constant with redshift. In addition, <J g is a small num- 
ber with a proportionally large uncertainty, and dividing 
by it contributes significantly to the scatter in V TO t/cr g . 
Nevertheless, due to the significant trends of ar g and V TO t 
with redshift, V ro t/o g shows strong evolution (4.8 o\ Fig- 
ure 6): 

log V rot /o g —1.55 = (— 0.48±0.10)(z-0.82)— (1.13±0.04). 

( 6 ) 

Tb determine whether the V rot of a galaxy influences 
our measurements of o g , in Figure 7 we examine how (J g 
evolves with redshift for galaxies from the DEEP2 Survey 
which are face-on in Hubble images (i.e., i < 30°) and in 
the same stellar mass range as the mass-complete sample. 
We find same evolution to within uncertainties for face- 
on galaxies: log <r g - 0,42 = (0.24 ± 0.04) (2 — 0.82) + 
(1.20±0.02), demonstrating that V roi does not influence 
our measurements of <r g . 

6. DOWNSIZING 

In this section, we demonstrate that much of the scat- 
ter in the relations of (t s /Sq . 5 and V Ta t/So.$ with red- 
shift is attributable to galaxy properties such as M, and 
size. We find that trends with M, and size are of similar 
strengths. 

First we consider M,. In Figure 8, as in Figure 6, we 
show o g /So.s and V rot /So .5 as a function of redshift for 
the mass-complete sample. However, in Figure 8, galax- 
ies axe divided into two M, bins: 9.8 < IogM»(Mg) < 
10.3 and 10.3 < logM»(M 0 ) < 10.7. Similar trends to 
those found for the entire mass-complete sample are also 
found separately for the low and high mass galaxies. It is 
also clear from Figure 8 that, at all redshifts, higher mass 
galaxies are more kinematically settled (i.e., higher V rat 
and lower <r 9 ) than lower mass galaxies. This is consis- 
tent with downsizing trends (e.g., Cowie et al. 1996), and 
we refer to this phenomenon as “kinematic downsizing.” 

Next we examine trends with size. We consider two 
size measurements: Hubble continuum sizes and emis- 
sion line extents, as described in §3. Hubble sizes are un- 
affected by seeing, and emission line extents are affected 
by seeing and are not corrected for it. Figure 9 is sim- 
ilar to Figure 8, except galaxies are color-coded accord- 
ing to size instead of M,. For both size measurements, 
three size bins are shown: < 3.5 kpc, 3.5-5 kpc, and > 5 
kpc. Figure 9 shows that on average, at all redshifts, the 
largest galaxies are the most kinematically well-ordered, 
the smallest the most disordered, and intermediate-sized 
galaxies have properties in between the large and small 
galaxies. This is the case for both size measurements. 
It demonstrates that the emission lines, which we use 
to measure galaxy kinematics, trace the galaxy continua 
well, and that we are able to measure reliable rotation 
velocities for the small galaxies in our sample. 

Figures 8 and 9 also demonstrate that the evolution in 
Figures 5 and 6 is not due to a special population of high- 
redshift galaxies with high values of o g and low values 
of Vrat ■ It occurs even if only massive/large galaxies are 
considered. 
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in Figure 6 , the quantities <?g / Ub.s end K-nt/ 1 S 0 5 are shown as a function of redshift;, but here galaxies are divided into two 
b ‘r v i° g M .* < ( low mas8 > black squares) and log M, > 10.3 (high mass, red triangles). Both mass ranges show similar trends 
with redshift as the full mass-complete sample, but the higher mass galaxies are the most well-ordered at all redshifts (high V rot /Sn s and 
low Og/So s)- Error bars and fits are the same as in Figure 5. 


7. SETTLING FRACTION 

lb better quantify how blue galaxies evolve in o g and 
Vrot, we use these parameters to define a kinematically 
settled disk galaxy, and study the fraction of settled disk 
galaxies as a function of M, and redshift. For the def- 
inition, we first look to the distribution of V TOi j<j g in 
Figure 3 and then to Hubble morphologies. Since we 
regard higher mass blue galaxies to be on average kine- 
matically settled disks at low redshift, we look to their 
values of V rat / cr g in Figure 3. A typical value for galaxies 
with lpg M.(Mq) > 9.8 is ~ 3 (i.e., log V rot /a g = 0.48). 
V e now consider Hubble morphologies. Kassin et al. in 
preparation show that Hubble morphologies are corre- 
lated With V ro t/(7 g such that galaxies with normal mor- 
phologies have higher values than galaxies with disturbed 
morphologies. We find that V r „t/a g ~ 3 provides a 
reasonable division between normal and disturbed mor- 
phologies, as demonstrated by the Hubble images of 6 
example galaxies from the K07 sample in Figure 10. 

Now that we have a quantitative definition of a set- 
tled disk galaxy (i.e., V rot /(T B > 3), we use it to study 
the fraction of blue galaxies which are settled, /settle- 
The evolution of f se ttu with redshift for the full K07 
sample divided into in four M, bins is shown in Fig- 
ure 11. All bins except for the lowest mass bin show 
decreasing f se ttie with increasing redshift (i.e., increas- 
ing f settle with time). It decreases from 38 — 24%, 
57 — 3l%, and 89 — 33% since z = 1.2 for mass bins 
of 3.0 < logAf,(M®) < 9.8, 9.8 < logM„(M Q ) < 10.3, 
and 10.3 < log M,(M Q ) < 10.7. These results are lower 
limits to the intrinsic evolution due to biases towards 
low/high mass galaxies in low/high redshift bins (Fig- 
ure 4). Furthermore, at all redshifts, the higher the M. 
of the population, the higher its f ae ttie - The same qual- 


itative behavior is found for settled galaxies defined as 
having values of V rot /cr g between 1 and 5. 

8 . KINEMATIC MEASUREMENTS IN THE LITERATURE 

In this section, we compare our measurements with 
those in the literature for local galaxies, and for galaxies 
at similar redshifts. We also compile kinematic measure- 
ments from the literature for galaxies at higher redshifts. 
Measurements of So. 5, V ro t, and o g from the literature are 
shown in Figure 12. AIL are from IFS observations. In the 
panels on the left and right, galaxies are color-coded ac- 
cording to stellar mass (for a Chabrier 2003 initial mass 
function) and emission-line extent, respectively, unless 
noted in the text below. Here we detail the measure- 
ments adopted from the literature in order of increasing 
redshift, and note if they differ from measurements used 
in this paper. 

To compare our results with local galaxies, a large 
homogenous sample is needed, and the GHASP survey 
(Epinat et al. 2010) provides just this. In Figure 12, we 
show the median values of So. 5 , V rat , and a g for galaxies 
in GHASP. Error bare on the median points show the rms 
scatter. GHASP spans a larger range in stellar masses 
than our sample: 9.0 < log < 11.7 for a ‘diet’ 

Salpeter IMF (Epinat et al. 2008), which is roughly 0.15 
dex higher than the Chabrier (2003) IMF used in this par 
per. Regardless, the median values are consistent with 
extrapolations of our relations to 2 = 0. 

In Figure 12, kinematics measurements for 68 galaxies 
from the IMAGES IFS Survey with a median redshift 
of 0.6 are shown (Fuech private comunieation, Puech 
et al. 2008, Neichel et al. 2009). The sizes used are 
continuum half-light radii. Although galaxies in IM- 
AGES are selected to be representative of galaxies with 
Af. > 1.5 x 10 10 M® at 0.4 < z < 0.75 (Yang et al. 
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FlO. 9.; — This figure is the same as Figure 8, except galaxies are divided into three size bins: < 3.5 kpc (black squares), 3.5-5 kpc 
(grey circles), and > 5 kpc (red triangles). Two sizes are considered: Hubble half-light radius (top; which is not affected by seeing), and 
emission line extent (bottom; which is affected by seeing and not corrected for its influence). Galaxies in size bins show similar trends 
to the full mass-complete sample. At all redshifts, the largest galaxies are the most well-ordered (i.e., high V ro t/So.6, low <r g /So.s), the 
intermediate-sized galaxies the next well-ordered, and the smallest galaxies the least well-ordered. Error bars and fits are the same as in 
Figure 5. For the Hubble continuum sizes, there is a deviation from the trends in the z — 0.7 bin because this is where we switch from 
using the F-band to the /-band for the size measurements. 


2098), similar to our selection, the median V TO t and o g 
of IMAGES galaxies are larger than ours at z = 0.6 by 
0.22 and 0.19 dex, respectively. It is possible that the 
median V rot is higher likely because the limited spatial 
coverage of the IFS may have caused V Tot to be over- 
estimated (Puech et al. 2008) . This especially seems to 
be the case for the smaller galaxies (< 3 kpc). The me- 
dian a g is higher because, unlike our procedure of fitting 
for V r ot and a g simultaneously, they are fit separately in 


IMAGES. Therefore, to derive an average a g from the 
velocity dispersion field of a galaxy, pixels affected by 
the seeing-induced peak at the center of a galaxy are 
ignored and a g is taken as the signal-to-noise weighted 
mean of the remaining pixels in the outer parts (Puech 
et al. 2007). This may result in a g measurements which 
are biased high because pixels with the highest signal-to- 
noise are ignored (e.g., Davies et al. 2011). 
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Fig. 10. — • Hubble/ ACS V + / images of settled (top) and unsettled (bottom) galaxies which are defined as having V TO t/a g > 3 and 
Vrct/dg < 3, respectively. The settled galaxies have normal morphologies, and the unsettled galaxies have disturbed morphologies (Kassin 
et al. in preparation). All galaxies are at z = 1 and have stellar masses in the range 10.4 < log < 11.1. The images are 6" on a 

side. 


8.1. Kinematic Measurements at 1.2 < z < 3.8 

Kinematics are difficult to obtain for galaxies at z > 
1.2 since emission lines are redshifted to near-infrared 
wavelengths where the sky is bright in continuum and 
line emission, making spectroscopy difficult. Neverthe- 
less, pioneering studies of galaxy kinematics at z > 1.2 
with single-slit spectrographs gave the first glimpse of the 
kinematic state of galaxies at these earlier times, albeit 
only the most highly star-forming systems (e.g., Erb et 
al. 2004, 2006). More recently, there has been a boom 
in IFS observations of galaxy kinematics at these red- 
shifts. Unfortunately, to obtain enough signal-to-noise 
to measure kinematics in a reasonable amount of observ- 
ing time, all samples at these redshifts are still biased 
toward galaxies with the brightest emission lines, which 
are the most highly star-forming systems. 

Many galaxies at z > 1.2 have kinematics which were 
unable to be fit by a rotating disk model (either due to. 
low signal-to-noise, small spatial extent, or a non-ordered 
velocity field). For these galaxies, we adopt as a rotation 
velocity in Figure 12 half of the velocity shear across the 
face of the galaxy ( | x V 3 h eaT ) . The quantity Vi/, Mr is the 
diiference between the minimum and maximum rotation 
velocities, and is not corrected for inclination. These 
galaxies are shown as small filled points in Figure 12. 

In Figure 12, kinematic measurements for 46 galaxies 


with a median redshift of 1.2 from the MASSIV Sur- 
vey are shown (Contini et al. 2012; Vergani et al. 2012). 
Galaxies in MASSIV were selected to have very bright 
emission lines, and a lower- limi t on emission linewidth 
was imposed (25A for 42 galaxies and 40A for 22 galax- 
ies; Contini et al. 2012). MASSIV galaxies have a median 
value of <r s which is consistent with our data, but higher 
median values of V rot and So.s- The difference is likely 
due to the sample selection. 

For 8 galaxies at z ~ 1.5, we adopt measurements 
from Wright et al. private communication for data in 
Wright et al. (2007) and Wright et al. (2009). All dual- 
component galaxies are treated as individual systems. 
Four of these galaxies were able to be modeled as rotating 
disks, and four were not. For galaxies which were able to 
be modeled, values of V TOt are not corrected for galaxy 
inclinations. For an additional 8 galaxies at z ~ 1.5, 
we adopt data from Lemoine-Busserolle & Lamareille 
(2010); all are modeled as rotating disks. For galaxies 
at z ~ 2, data are adopted from the SINS Survey via 
private communication from N. Forster-Schreiber for the 
most recent models and data for galaxies in Genzel et 
al. (2008), Forster-Schreiber et al. (2009), and Cresci et 
al. (2009). A total of 14 SINS galaxies were able to be 
modeled as rotating disks, and 33 were not. For galaxies 
from SINS which were not able to be modeled, we esti- 
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Pig. 11.— The fraction of galaxies in the full K07 sample which 
are settled (defined as V rot /<7g > 3), f, e ttu, is shown as a func- 
tion of redshift for different stellar mass bins. The quantity f settle 
decreases with redshift (increases with time) for all but the low- 
est mass bin. At all redshifts, the higher the stellar mass of the 
papulation, the higher its f settle- 


mate o g as the line-width less the velocity shear (i.e., v 0 be 
from Forster-Schreiber et al. 2009 - V S hear)- For an addi- 
tional 16 systems at z ~ 2, none of which were able to be 
modeled as a rotating disk, we adopt data from Law et 
al. (2009). Three of the galaxies from Law et al. (2009) 
which have dual-components are considered as separate 
systems here. For 4 galaxies at z ~ 3, data are adopted 
from Lemoine-Busserolle et al. (2010); 3 are modeled as 
rotating disks and 1 is not. In addition, at z ~ 3, 12 
gflaxies which were modeled as disks and 7 which were 
not, are adopted from from the AMAZE/LSD Survey via 
private communication from A. Gnerucci for data from 
Gnerucci (2011) and Troncoso et al. in preparation. 

Summary 

For measurements in the literature over 1.2 < z < 3.8, 
there is a large scatter in kinematic parameters and no 
clear trends with redshift. This is likely due to sample 
selection since currently only the brightest systems at 
these redshifts can be observed in reasonable amounts of 
observing time. In addition, many of the higher redshift 
galaxies are expected to merge and become early-type 
galaxies by at least today, inducing a “progenitor bias.” 

These data demonstrate that there exist galaxies at 
earlier times with rotation velocities similar to local 
galaxies. However, they differ from local galaxies due 
to large values of a g (e.g., Forster-Schreiber et al. 
2009; Lemoine-Busserolle & Lamareillc 2010; Lemoine- 
Busserolle et al. 2010; Gnerucci 2011; Vergani et al. 
2012). The data from the literature also show the ex- 
istence of a population of higher redshift galaxies with 
small sizes and little or no discernible rotation, but with 
significant values of <r g and generally small sizes (Law et 
al. 2009: Forster-Schreiber et al. 2009). Although some 
observations at z ~ 2 may be too shallow to detect the 


full rotation gradients (Law et al. 2009; Forster-Schreiber 
et al. 2009), until these smaller galaxies earn be observed 
with significantly longer exposure times to detect pos- 
sible faint disks, whether they have a significant V TOt is 
unknown. Of the galaxies modeled as rotating disks at 
1.2 < z < 3.0, the more massive systems generally have 
higher values of V rot and 5o.s, but no trends with c g . 
At z a 3, nearly all galaxies are small, but have large 
rotation velocities. 

9. CONCLUSIONS 

We use a sample of 544 galaxies to study the inter- 
nal kinematics of blue galaxies over the last ~ 8 billion 
years, and find significant evolution. Galaxies become 
progressively more well-ordered with time as disturbed 
motions decrease and rotation velocities increase. In ad- 
dition, galaxy potential well depths increase -with time. 
At all redshifts, the more massive galaxies are on average 
the most kinematically evolved. The lowest mass galax- 
ies in our survey, which are observed at low redshift, are 
not well-ordered. They may be settling today, or may 
never become settled disks. These trends with mass are 
consistent with downsizing trends for other galaxy prop- 
erties, and we refer to this phenomenon as “kinematic 
downsizing.” 

We define a kinematically settled disk galaxy as hav- 
ing a ratio of ordered to disordered motions ( V ra t/<7 g ) 
of 3. We find that the fraction of settled disk galaxies, 
/settle, increases with time since z — 1.2 for all but the 
lowest mass bin in our survey, which we can only probe 
at low redshifts. The fraction f settle is larger the more 
massive the galaxy population is, and this is independent 
of redshift. These results do not change qualitatively if 
a settled disk galaxy is defined with values of V TO t/a g 
between 1 and 5. 

We speculate that this evolution may be due to a de- 
creasing rate of minor mergers/mass accretion which is 
expected in a ACDM Universe (e.g., as invoked by Oser 
et al. 2012, to explain the evolution of early-type galax- 
ies) and may be linked to the precipitous decrease in 
galaxy star-formation rate and molecular gas fraction 
since z = 1.2. Phenomena which are the causes and ef- 
fects of star-formation, such as minor mergers/accretion 
and gas motions, may increase disturbed motions. In ad- 
dition, the large molecular gas reservoirs of higher red- 
shift galaxies (Daddi et al. 2010; Tacconi et al. 2010) 
provide fuel not only for for star-formation, but also 
for disturbed motions through the Toomre Q parame- 
ter (Toomre 1964) as in e.g., Silk & Norman (2009). In 
addition, kinematic downsizing may be an effect of star- 
formation which began earlier and declined more quickly 
in massive galaxies, compared with a later on-set of star- 
formation in less massive systems. And this might be 
related to the minor merger rates for high and low mass 
galaxies. 

We compile measurements of galaxy kinematics from 
the literature over 0 < z < 3.8. Observations of the 
kinematics of local disk galaxies reveal few examples of 
massive disk galaxies where disturbed motions play as 
an important of a role as in similar mass disk galaxies 
at higher redshifts (e.g., Epinat et al. 2010). Our com- 
pilation of measurements over 1.5 < z < 3.8 shows no 
trends with redshift, likely due to sample selection, since 
only the most highly star-forming systems at these red- 
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Fig. 12. — These plots are the same as in Figure 5, except here the z axes are expanded to include data from the literature for lower and 
higher redshift galaxies. In the plots on the left and right, galaxies are color-coded according to M '« and emission line extent, respectively. 
For galaxies the local GHASP survey, only the medians are shown (open black triangles), with error bars depicting the rms scatter. 
Otherwise, large open symbols denote individual galaxies where disk models were fit. Small filled points denote individual galaxies which 
were unable to be fit by a rotating disk model; for these galaxies, V Bhear and an average tr g arc plotted. Large filled symbols represent 
survey medians. The providence of each galaxy from the literature is indicated by the symbol shape (regardless of whether it is open or 
filled), and references are labeled in the central panel on the left. 
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shifts can currently be studied. However, it is clear that 
most of these higher redshiffc systems have a significant 
amount of disturbed motions, whether or not they show 
evidence of rotation. 

The exact causes of the disordered motions (as quan- 
tified by o g ) are not yet known. They may indicate a 
formation mechanism for the thick disk, or a gas disk 
that is thickened by stirring, gas motions due to star- 
formation, and/or minor mergers. Perhaps we are wit- 
nessing the phenomenon of disk settling as disturbed mo- 
tions are turned into ordered rotation via mechanisms 
such as circularization of gas orbits and/or gaseous dis- 
sipation to the disk plane. We do not yet know how this 
process of settling to a disk occurs, and whether bary- 
onic dissipation is more similar to a smooth process or 
a repeated minor merging of clumps. Only with very 
high spatial resolution kinematic measurements, as can 


be obtained with IFS and adaptive optics on future 30- 
40-meter telescopes or the James Webb Space Telescope, 
will we be able to unambiguously reveal the nature of 
the processes which cause o g . 
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